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Abstract 

This work results from our attempts to solve Boltzmann-Sinai's 
hypothesis about the ergodicity of hard ball gases. A crucial element 
in the studies of the dynamics of hard balls is the analysis of special 
hypersurfaces in the phase space consisting of degenerate trajectories 
(which lack complete hyperbolicity). We prove that if a flow-invariant 
hypersurface J in the phase space of a semi-dispersing billiard has a 
negative Lyapunov function, then the volume of the forward image 
of J grows at least linearly in time. Our proof is independent of the 
solution of the Boltzmann-Sinai hypothesis, and we provide a complete 
and self-contained argument here. 

Keywords: Hard balls, Bolzmann-Sinai hypothesis, semi-dispersing billiards, 
ergodicity, Lyapunov function. 

1 Introduction 

The ergodic hypothesis for gases of hard balls was first put forward (in rather 
vague terms) by L. Boltzmann back in the 1880 's, and then formalized by 
Ya. G. Sinai in 1963 [Sin63 . It states that the gas of > 2 identical 
hard balls (of small radius) on a torus Tor'^, > 2, is ergodic, after certain 
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necessary reductions. The latter mean that one fixes the total energy, sets the 
total momentum to zero, and restricts the center of mass to a certain discrete 
lattice within the torus. The assumption of a small radius is necessary to 
have the configuration space connected. 

The motion of N hard balls in a d-dimensional torus naturally reduces 
to that of a single billiard particle in a {dN — (i)-dimensional domain P 
with specular refiections off its boundary dT>. The domain "D, which is the 
configuration space of the system, is obtained by the removal of A^(A^ — l)/2 
cylinders from the torus Tor°'^~'^ (each cylinder corresponds to a pair of 
colliding balls). Thus the boundary &D is convex (as seen from inside the 
domain V), but not strictly convex, which makes the billiard in V semi- 
dispersing. This is true, unless = 2, in which case we have a single 
'cylinder' that actually becomes a sphere, thus dT> is strictly convex, which 
makes the corresponding billiard dispersing. 

The system of iV = 2 balls (disks) in dimension d = 2 was thoroughly 
investigated by Sinai |Sin70j : he proved its hyperbolicity, ergodicity, and 
K-mixing. In this seminal paper Sinai also developed a general theory of 
planar dispersing billiards. Then Sinai and Chernov |SC87j extended these 
results to systems of = 2 balls in any dimension c? > 2, as well as to other 
multidimensional dispersing billiards (there was a notable oversig htininniTi 
that was corrected later in |BCST02] ). 

The dynamics of hard balls (and related billiards) is not entirely smooth 
- singularities occur at 'grazing' collisions of two balls and at multiple col- 
lisions of > 3 balls (if they hit each other simultaneously). These events 
induce hypersurfaces in the phase space consisting of singular trajectories; 
they are called singularity manifolds and cause major troubles in the analy- 
sis. In addition, there is an issue of uniform versus nonuniform hyperbolicity. 
Dispersing billiards are characterized by strong hyperbolicity with uniform 
expansion and contraction rates. On the contrary, semi- dispersing billiards 
are only non-uniformly hyperbolic, so that their expansion and contraction 
rates may be arbitrarily slow. This constitutes a crucial difference between 
systems of = 2 (dispersing) and N > 2 (semidispersing) cases, the latter 
are much harder to deal with. 

The first major result for systems of > 2 balls was the theorem on 
'local' ergodicity by Sinai and Chernov |SC87j (a plan of its proof was earlier 
laid out by Sinai |Sin79j ): loosely speaking it states that for any hyperbolic 
phase point x there is an open neighborhood U{x) that belongs to one ergodic 
component (mod 0). That theorem was proven under an assumption later 
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referred to as Chernov-Sinai Ansatz |KSS90t IKSSQlj : roughly speaking it 
states that typical points on singularity manifolds must be hyperbolic. 

A strategy for proving complete hyperbolicity and 'global' ergodicity for 
systems of > 2 hard balls was first proposed by Sinai and Chernov in an 
unpublished manuscript |SC85j and further developed and refined by Kramli, 
Simanyi, and Szasz |KSS90t IKSS9H IKSS92j . It is based on the following 
observations. 

A phase point x fails to be hyperbolic in two possible ways. First, it may 
happen that along the trajectory of x the system of A^ balls splits into two (or 
more) groups (clusters) of balls that only interact within clusters. As unlikely 
as it seems, such an anomaly does occur, it is called splitting. Second, even 
without splitting the balls may collide in such a 'degenerate' manner that 
expansion and contraction of tangent vectors only occur in a tangent subspace 
of some lower dimension. The simplest example is obtained by letting all the 
A^ balls move along a closed geodesic in the torus (with different but parallel 
velocity vectors). In such examples the balls seem to 'conspire' to prevent 
complete hyperbolicity; in any case their positions and velocity vectors must 
satisfy very stringent requirements. We call this phenomenon degeneracy. 
Thus the complete hyperbolicity can only be destroyed by either splitting or 
degeneracy. 

In order to prove the hyperbolicity for the system of A^ balls, it is enough 
to verify that splitting and degeneracy occur on subsets of zero measure. 
To establish ergodicity, one needs more than that - the set of 'bad' phase 
points (where splitting or degeneracy occur) must not separate two ergodic 
components. For instance, it is enough to show that the set of all bad phase 
points has topological codimension two or more. In addition, one needs to 
verify the Ansatz to utilize the local ergodic theorem. All these proofs are 
carried out inductively, by assuming the ergodicity of systems of n hard balls 
for all n < N. 

The above strategy was successfully employed by Kramli, Simanyi and 
Szasz who proved hyperbolicity and ergodicity for A^ = 3 balls in any di- 
mension |KSS91j and for A^ = 4 balls in dimension d > 3, see |KSS92j . 
Then Simanyi jSim92a| ISim92bj proved hyperbolicity and ergodicity when- 
ever N < d (this covers systems with an arbitrary number of balls, but only 
in spaces of high enough dimension, which is a very restrictive condition). By 
using the induction hypothesis mentioned above, Simanyi |Sim92al ISim92bj 
also found a general argument showing that in any system of hard balls, split- 
ting only occurred on a subset of zero measure and topological codimension 
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two (or more) in phase space, thus sphtting was no longer a problem. But 
degeneracies remained a key problem - in all the above papers they were 
handled by a direct analysis involving various cases specific for every system 
of hard balls. This 'case study' became overly complicated with every extra 
ball added to the system, and it was clearly impractical for > 4 balls. 



Further progress required novel ideas, and Simanyi and Szasz |SS99j em- 
ployed the methods of algebraic geometry. They assumed that the balls had 
arbitrary masses mi, . . . ,mAr (but the same radius r). Now by taking the 
limit mjv — i> they were able to reduce the dynamics of balls to the mo- 
tion of — 1 balls, thus utilizing a natural induction on N. Then, algebraic 
methods allowed them to effectively analyze all possible degeneracies, but 
only for typical (generic) vectors of "external" parameters (mi, . . . ,m7v,T); 
the latter needed to avoid some exceptional submanifolds of codimension 
one in R^"*"^, which remained largely unknown. This approach led to a proof 
of complete hyperbolicity (but not yet ergodicity) for all iV > 2 and 



d > 2, and for generic (mi, . . . ,m]\f,r). Later Simanyi alone upgraded the 
hyperbolicity to ergodicity for all systems of > 2 hard balls and typical 
(generic) values of external parameters (mi, . . . , m^, r); this was done in two 
separate papers: |Sim03j covered the case d = 2 and |Sim04ij dealt with d > 2. 
Both papers again use algebraic methods, thus their results are restricted to 
'generic' masses. This was not quite satisfactory as the most interesting case 
of equal masses mi = ■ ■ ■ = (it is exactly this case that motivated both 
Boltzmann and Sinai) remained open, as there was no guarantee that the 
exceptional submanifolds would avoid it. 

It is clear that methods of algebraic geometry, however powerful, have to 
be abandoned if we are to obtain hyperbolicity and ergodicity for all gases 
of hard balls (including the case of equal masses). Recently Simanyi |Sim02j 
made a partial progress in this direction: he employed purely dynamical 
arguments to establish complete hyperbolicity (but not yet ergodicity) for all 
systems of hard balls (including equal masses mi = ■ ■ ■ = 171^). In the course 
of that work, Simanyi showed that degeneracies occurred on a countable 
union of smooth submanifolds in the phase space with codimension at least 
one. This implies, of course, that degeneracies are restricted to a subset of 
zero measure, which is enough for hyperbolicity, but ergodicity would require 
those submanifolds be of codimension at least two. 

In an attempt to find a general dynamical argument that would rule 
out the existence of degeneracies on submanifolds of codimension one (hy- 
persurfaces) we showed that if such hypersurfaces existed, then they would 
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necessarily grow in size (volume) either in the past or in the future (depend- 
ing whether the Lyapunov function on them is positive or negative). This 
fact turned out to be essential in a subsequent proof of ergodicity by one of 
us (Simanyi), under certain conditions. Here we prove the volume growth of 
hypersurfaces, as a separate fact, in the context of general semi-dispersing 
billiards (which of course include all systems of hard balls) . 



2 Results and proofs 

Let P be a bounded connected domain in R"' or in the torus Tor*^; the bound- 
ary dV is a finite union of smooth compact hypersurfaces (each, possibly, 
with boundary): 

91? = Ti u ■ ■ • u r^. 

A billiard system is generated by a pointwise particle; it moves freely (with 
a constant velocity vector) at unit speed in V and gets specularly reflected 
at the boundary dV (by the classical rule "the angle of incidence is equal to 
the angle of reflection" ) . 

The phase space of the billiard system is the unit tangent bundle VL = 
V X S'^~^ over T). Thus phase points are pairs x = {q, v) where g G I? is the 
position and v G S'^~^ the unit velocity vector of the billiard particle. The 
billiard dynamics generates a flow : Q Q. It is a Hamiltonian flow that 
preserves its Liouville measure fi; the latter is a direct product of uniform 
measures on V and S"^~^. 

At every (regular) boundary point g G Fj \ dTi, 1 < z < r, we denote by 
z/(g) the (unique) unit normal vector to Fj pointing into T>. The curvature 
operator (the second fundamental form) JCq is a self-adjoint linear transfor- 
mation acting on the tangent space TgYf, is given (to the linear order) by 

u{q + 6q) = iy{q) + ICq{6q). 

The billiard in T> is said to be dispersing (semi-dispersing) if JCg is positive 
definite (resp., positive semi-definite) at every regular point q G dV. Geomet- 
rically, this means that &D is strictly convex (resp., just convex) as seen from 
inside V. For a recent detailed studies of the dynamics in semi-dispersing 
billiards we refer the reader to |BCST03] . 

Now let J C I? be a small smooth compact hypersurface (a submanifold 
of CO dimension-one with boundary), which is locally flow- invariant, i.e. for 
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every x G J \ dJ there is e > such that $*x G J for all \t\ < e. For every 
X = {q,v) G J \ dJ we denote by = {z, w) G %Q a non-zero normal 
vector to J, i.e. such that for any tangent vector {6q, 6v) G %fl the relation 
{6q, Sv) G J is equivalent to {6q, z) + {6v, w) = 0; here ( ■ , ■ ) denotes the 
Euclidean inner product in W^. 

Observe that {w, v) = because ||f || = 1 for all x = (g, v) G Q. Also since 
J is flow-invariant, we have (f , 0) G 7^ J, hence {z,v) = 0. Thus w,z E v^, 
where denotes the hyperplane in Mf^ orthogonal to v. 

The value Q{nx)- = {z,w) is called the (infinitesimal) Lyapunov func- 
tion, see jKB94j or part A. 4 of the Appendix in |Ch94j . For a detailed ac- 
count of the relationship between the Lyapunov function Q, the symplectic 
geometry in Q, and the dynamics, see |LW95j . 

Remark. Since the normal vector = [z, w) to J is only determined up 
to a nonzero scalar multiple, the value Q{nx) is only determined up to a 
positive multiple. However, in our considerations the sign of Q{nx) will be 
most important, and this sign is clearly independent of the choice of n^- 

Next, for any t G M the image Jt = ^^(J) is also a locally flow-invariant 
compact hypersurface in Q. Given a point yo = {qo, Vq) E J \ dJ, its image 
Vt = {qt,vt) = $*(l/o) is a regular point in Jt (i.e. yt ^ dJt), unless yt is a 
reflection point. Given a normal vector uq to J a.t yo, we define the unique 
normal vector nt to Jt at yt by 

nt = 

i.e. 

(2.1) {Dyy{Sy),nt) = {6y,no) W 6y G Ty.Q 

It is enough to require this for all 6y = {6q,6v) G Ty^^ such that both 6q 
and Sy belong in Vq, because w,s E as well. 

We develop explicit formulas for rit (t > 0). If there is no collisions 
during an interval of time [s,t] on the orbit {yt}, then the relation between 
{6qs, 6vs) G Ty^Q and {6qt, 6vt) = Dy^^^'^lSqs, 6vs) is obviously 

6vt = 6vs, 

6qt = Sqs + (t- s)6vs, 

from which we obtain that for all 6ys = {Sqs, Svg) E Ty^Q 

{Sys,ns) = {5qt - (t - s)5vt, Zs) + {5vt, Ws) 
= {Sqt, Zs) + {5vt, Ws-it- s)zs). 
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hence 



(2.2) 



fit = {Zt, Wt) 



[Zs, Ws-{t- s)Zs). 



Observe that 



(2.3) 



Q{nt) = Q{ns) 



{t-s)\\zsf <Q{n,). 



Next let the orbit {yt} coUide (reflect) at a boundary point q G dV at 
some time t > 0. The velocity vector v~ = Vt-o is transformed to = Vt+o 
by the rule 



where (p is the angle between the outgoing velocity vector f + and the unit 
normal vector z/(g) to dT>, i.e. cos 9? = {v~^, i^iq))- 

At the collision, the flow transforms tangent vectors {Sq~ , 6v^) G Ty^_^^Vt 
(we can assume that both 5q~ and 5v~ belong in (f ")"*", see above), according 
to the standard rules 



where the operator R: TqT> TqT> is the orthogonal reflection across the 
tangent hyperplane Tg&D, V: (f")"*" TqdT> is the f "-parallel projection of 
iy^)'^ onto TqdV, V* is the adjoint of V, i.e. it is the z/(g)-parallel projection 
of TgdV onto (f")"*", and K. is the curvature operator of dV at q (see above). 

The formulas (|2.4|) are given, e.g., in |SC82t Section 1] or |KSS90t Propo- 
sition 2.3]. We can also rewrite them as 



where Vi is the f+-parallel projection of (f ^)"'" onto TqdT>. 

Now let n~ = {z~,w~) denote the normal vector to Jt at yt-o, i.e. right 
before the collision, then 



{{6q", 6v'),n~) = {R{6q+), z') + {R{bv^) ~ 2cos(f)RV*ICVi{6q+), w~) 



Thus the normal vector n+ = {w^.z'^) to Jt at y^+o (after the collision) is 



= V + 2 cosy? z/(g). 



(2.4) 




R{5q-), 

R{6v') + 2 cos ^RV*)CV{6q-), 



6q- = i?(V), 

Sv- = R{Sv+) - 2cos0 W;/CVi(5g+), 



R{z-) -2cos(f)V{}CViR{w-)) + {6v+, R{w-)). 



n 



+ 



{R{z-) -2cos4>V*ICViR{w-), R{w 



))■ 
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Observe that \\w~^\\ = \\w \\ and 

Q(n+) = Q{n-) -2 cos (j){V*JCViR{w"), R{w-)) 

= Q{n-) -2 cos (f){K:ViR{w-), ViR{w-)) < Q(n"), 

because /C > for semi-dispersing billiards. 
We summarize our results: 

Proposition 2.1. The Lyapunov function Q{nt) is a monotone non-increasing 
function oft. The norm \\wt\\ is a continuous function oft. The norm \\zt\\ 
is a piece-wise constant function oft changing only at collisions. 

Our normal vector Ut is related to the rate of expansion of the hypersur- 
face J by the flow $*: 

Proposition 2.2. The flow $* expands the volume of the hypersurface J at 
a point yQ E J by a factor ||^^t||/||?^o|| ■ 

Proof. According to (|2.H) . for every 5yt E Ty^Q we have 

{Syt, rit) = {Dy^<^-\6yt), uq) . 
Substituting 6yt = rit gives 

ll^tll _ {Dy,^'\nt), no) 
\\no\\ \\nt\\ \\no\\ 

A simple geometric inspection shows that the right-hand-side is exactly the 
factor of linear expansion by $"* between the points yt and yo in the direction 
transversal to Jt, and since preserves volume in Q, it is the factor of volume 
expansion of the hypersurface J by between the points yo and yt- The 
proposition is proved. □ 

From now on we assume that Q{no) < 0. This assumption is motivated 
by the following considerations (which will not affect our further arguments, 
though). If there is a flow- invariant compact hypersurface J G Q, it might 
divide the phase space Q into two (or more) non-interacting open subsets, 
thus preventing the ergodicity of the flow. One wants to show that such 
hypersurfaces cannot exist. There is a separate argument |Simn8| Remark 
7.9] showing that the Lyapunov function cannot stay identically zero on the 
entire trajectory of any non-splitting phase point, thus Q < or Q > on 
parts of J. We will show that if Q < 0, then the volume of J is expanded un- 
der the flow, at least linearly in time, hence J cannot be compact. Similarly, 
if Q > 0, then the volume of J is expanded in the past. 
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Proposition 2.3. If Q{nQ) < 0, then Q{nt) is a monotone decreasing func- 
tion oft, and the norm \\wt\\ of the 'velocity' component of the normal vector 
nt = {zt, Wt) is a monotone increasing function oft. 

Proof. We have Q{ns) = {zs, Ws) < Q{no) < for all s > 0. Hence ||zs|| 7^ 0, 
and so Q{nt) < Q{ns) for any collision-free time interval due to ()2.3|) . 
Also, by 

Ikif = Iksf - 2(t - s) Q{ns) + (t - s)^||2;,f > ||w,f , 

and \\wt\\ does not change at collisions. □ 

Proposition 2.4. IfQ{no) < 0, then the function \\wt\\/\Q{nt)\ is monotone 
decreasing in t. 

Proof. For any s > we have 

Q{ns) , j_ 
Ws = -r. — p Zs + W^, 

1 1 1 1 

where is a vector orthogonal to Zg, hence 



Ws 



1 \\w„ 



±112 



1/2 



Due to ()2.2j) . for any collision- free time interval {s,t) we have Zt = Zg and 
= w^, while |Q(ns)| < \Q{nt)\. At colhsions, \\wt\\ stays constant, while 
|(5(^t)| instantaneously increases, thus the required property holds. □ 

Proposition 2.5. If Q{no) < 0, then for all t >0 

II ||.|| II , \Qino)\t 

\\wt\\ > \\wo\\ H —. 

IFoII 

Proof. For any t > 0, which is not a collision time. 



dt dt \\wo\\ 

where we used fl2.2|) and the previous proposition, hence 

dt ~ \\wo\\ 

Since \\wt\\ does not change at collisions, the proposition follows. □ 
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Next we assume that Q{t) < on all of the manifold J; more specifically 
(2.5) Q{no) = {zo, wq) < -cq < 

holds true uniformly across J for a unit normal vector field no{y) {y G J) 
with a fixed constant cq > 0. 

Theorem 2.6. // holds, then $* expands the volume of the hypersurface 
J at any point y E J by a factor > 1 + Cgt for all t > I/cq. 

Proof. Combining Propositions 12.21 and 12.51 gives 

- Cot 

At = \\nt\\ > \\wt\\ > \\wo\\ + - — r-. 

IFoII 

Since \\wo\\ < WuqW = 1, this lower bound is at least 1 + cot for t > I/cq. □ 
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